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ABSTRACT 

Four close radio sources in the International Celestial Reference Frame (ICRF) catalog were 
observed using phase referencing with the VLBA at 43, 23 and 8.6 GHz, and with VERA at 23 
GHz over a one year period. The goal was to determine the stability of the radio cores, and to 
assess structure effects associated with positions in the ICRF. Although the four sources were 
compact at 8.6 GHz, the VLBA images at 43 GHz with 0.3-mas resolution showed that all were 
composed of several components. A component in each source was identified as the radio core 
using some or all of the following emission properties: compactness, spectral index, location at 
the end of the extended emission region, and stationary in the sky. Over the observing period, the 
relative positions between the four radio cores were constant to 0.02 mas, the phase referencing 
positional accuracy obtained at 23 and 43 GHz among the sources, suggesting that once a radio 
core is identified, it remains stationary in the sky to this accuracy. Other radio components in 
two of the four sources had detectable motion in the radio jet direction. Comparison of the 23 
and 43 GHz VLBA images with the VLBA 8.6 GHz images and the ICRF positions suggests that 
some ICRF positions are dominated by a moving jet component; hence, they can be displaced 
up to 0.5 mas from the radio core, and may also reflect the motion of the jet component. Future 
astrometric efforts to determine a more accurate quasar reference frame at 23 and 43 GHz and 
from the VLBI2010 project are discussed, and supporting VLBA or EVN observations of ICRF 
sources at 43 GHz are recommended in order to determine the internal structure of the sources. 
A future collaboration between the radio (ICRF) and the optical frame of Gaia is discussed. 

Subject headings: astrometry, quasars: general, radio continuum: galaxies, surveys 

1. Motivation for the Experiment errors for many of these sources is less than 0.06 

mas obtained by averaging decades of observa- 
The Second Realization of the International Ce- ^^^^^^ ^nd they define the axes of a radio-based 

lestial Reference Frame (ICRF2) is defined by the quasi-inertial frame to about 0.01 mas precision, 

positions of 295 compact radio sources located ^his positional accuracy is presently limited by a 

around the sky (jFev et al 1 120101 ) . The positional 
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typical residual variable tropospheric refraction of 
about 20 psec ( 0.15 mas) for a time-scale of tens of 
minutes and angular size scale of tens of degrees. 

The radio source emission for most sources 
in the ICRF2 catalog (quasars and galaxies) is 
dominated by a strong compact radio component 
(called the radio core) with a size < 0.1 mas. This 
radio core is probably located at the base of a ra- 
dio jet that is formed perpendicular to the accre- 
tion disk about 1000 to 10000 Schwarzschild radii 
from the m assive object in the galactic nucleus 
( Marscherl 12006.) . For quasars and radio galaxies 
at cosmological distances, this offset is less than 
about 0.01 mas, so the radio core should represent 
a fixed point in the sky at this precision. However, 
there is often more extended emission (generically 
called the jet) that may be connected to the core, 
or consist of discrete components. The emission of 
the core and jet both vary with time, and the jet 
often contains components that move away from 
the core. Because the emission near the core be- 
comes more opaque at lower frequencies, the de- 
termination of the core position becomes problem- 
atic. Although the magnitude of this uncertainty 
caused by structure is less than that from tropo- 
sphere refraction, it is becoming an important part 
of the position error budget as the tropospheric 
modeling and measurements improve. 

In order to determine the effect of structure and 
its changes on the apparent radio core positions 
obtained from radio observations, the effect of the 
tropospheric refraction variations must be signif- 
icantly decreased in order to reach the structure 
astrometric effects. By using a VLBI technique 
called phase referencing, where several sources 
that are closely spaced on the sky are observed 
simultaneously or in quick succession, a large part 
of the tropospheric refraction can be canceled. For 
frequencies as high as 43 GHz, the needed phase 
stability (residual tropospheric-induced and other 
phase errors should be less than about 1 radian) 
can be obtained most of the time if the separations 
among the sources are less than a few degrees, and 
the switching times between sources is less than 
about one minute. 



2. Observations and Initial Data Process- 
ing 

2.1. Source Selection 

The source selection was made to accommodate 
observations by the VLBI Exploration of Radio 
Astrom etry (VERA) four-e lement VLBI array in 
Japan ( Honma et al I [ioosh that can observe two 
sources simultaneously if they are separated by no 
more than 2.2°. A search of the 8.6-GHz ICRF 
catalog found several sets of sources and the best 
set for phase referencing were the four sources: 
0547-^234, 0554-^242, 0556-f-238 and 0601+245. 
The relative position of the four sources are illus- 
trated in Fig. [TJ The positions and approximate 
intensity parameters for the sources are shown in 
Table [1] The apriori positions, taken from the 
ICRF catalog, are fiducial locations for each source 
(phase center used in the correlation of the data) 
to which the measured offset positions and struc- 
tures are referred. The approximate correlated 
fiux densities at each frequency are given for the 
projected spacing of 200 km (S) and at 5000 km 
(L). The values have been taken from these obser- 
vations. All of the sources are somewhat variable. 
Because the sources are within 5° of the galactic 
plane, there are no firm identifications, and their 
redshifts are not known. The sources are prob- 
ably quasars with a redsh ift between 0.5 to 2.0 



(iHardwick fc Schade lIToOOl) 



The source 0556-1-238 is one of the 295 defining 
sources in the ICRF2 catalog. The other three 
sources are in the extended ICRF2 list, but are 
not one of the defining sources because of their 
somewhat limited use in the ICRF sessions over 
the past thirty years. 

Although this quartet of sources was chosen be- 
cause of their proximity in the sky and their rel- 
ative strength and compact structure at 8.6 GHz, 
other deeper properties of these sources — internal 
structure and its frequency dependence, and vari- 
ability time-scale — were unknown at the higher 
frequencies and resolutions used for these obser- 
vations. These four sources are a representative, 
although limited, sample of good quality ICRF 
sources for the deeper investigations reported in 
this paper. 
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2.2. Observing Sequences 

The list of observing sessions are given in Ta- 
ble [5] The initial observing request was submitted 
to the Japanese VERA array, and the first ob- 
serving session at 23 GHz occurred on 18 April 
2008. The session was 6.5 hours in length, and 
cycled among the source pairs: 0547-0554; 0547- 
0556; 0554-0556; 0554-0601; 0556-0601, ah with a 
separation less than 2.2° so VERA could observe 
each source-pair simultaneously. Only the source 
pair separation 0547-1-234 to 0601-1-245 was over 
this angular limit. Each pair was observed for an 
8-minute scan, and then each of the five pairs had 
8 such observations over each session. These data 
enabled accurate images and pair separations to 
be determined. After the second VERA session on 
24 May 2008, an observing proposal was submit- 
ted to the VLBA in order to observe at somewhat 
higher resolution and at 43 GHz since the prelim- 
inary VERA results indicated the complex nature 
of the sources. 

Six VLBA sessions (see Table [2] at 23 GHz and 
43 GHz were scheduled between 20 Dec 2008 and 
12 Dec 2009. For the most accurate phase refer- 
encing stability, a primary reference source that is 
observed every other scan was chosen. The source 
0556-1-238 was chosen as the primary calibrator 
because it was the most compact (ratio of L/S 
from Table [T] was largest), and it was relatively 
strong. This permitted the monitoring of fast 
phase changes caused by the tropospheric refrac- 
tion variations. The other sources were then al- 
ternated observed in turn with the primary phase 
calibrator. This observing/reduction scheme uses 
the primary reference source to determine the in- 
strumental/atmospheric/astrometric errors in its 
direction and assumes that the same errors apply 
to the other sources (targets). The positional ref- 
erence frame is also tied to the apriori position of 
the primary reference source. 

The observing sequence used was 0556-0547- 
0556-0554-0556-0601-0556-0547...). Each scan 
was 18 sec in length, with 10 sec needed for source 
switching. The separation between the center of 
two 0556-1-238 scans was only 55 sec, sufficiently 
close in time to maintain phase coherence at 43 
GHz except under adverse weather conditions or 
at very low elevations. The basic sequence took 
3 min, and ten sequences were made over a 30 



minute period. These 30-minute blocks were then 
alternated between 23 GHz and 43 GHz observing 
frequencies. 

For 40 minutes near the middle of each VLBA 
session, 1-min scans at 23 GHz of strong ICRF2 
sources over the sky were made in order to deter- 
mine the residual tropospheric refraction model 
in addition to that used for the baseline corre- 
lator model. This observing technique is calle d 



"DELZN", and is described bv lReid et aTl (|2009f ) 



The final VLBA session on 31 Jan 2010 for one 
hour at 8.4 GHz, used the same scan cycle in order 
to obtain the relative source positions and struc- 
ture at this lower frequency at the end of the mon- 
itoring period. 

2.3. Data Correlation and A-priori Pro- 
cessing 

The correlator models used for the VERA and 
VLBA observations used accurate locations for the 
antennas, up-to-date earth orientation (EOF) pa- 
rameters, and initial tropospheric delay models. 
The positions listed in Table [T] are the assumed 
positions (phase center) for the sources for all ses- 
sions. These positions are merely the reference 
point from which offset positions are determined. 
They are often close to the location of the maxi- 
mum intensity of the radio source. 

For all VERA sessions, the correlated data were 
sampled every 2 sec, with 512 channels in one 
frequency band, covering at total bandwidth of 
128 MHz, centered at 23.564 GHz. After flag- 
ging data associated with mal-functioning anten- 
nas or system temperatures in excess of 200K (the 
nominal operating system temperature of 100- 
120K increased significantly during periods of poor 
weather), the remaining data were then averaged 
to 6-sec intervals and 32 channels covering the 
128 MHz bandwidth. No bandpass corrections in 
phase and amplitude were needed. Each of the 
five source pairs for each session were processed 
separately: 0547-0556(PR), 0554-0556(PR) 0601- 
0556(PR), 0547-0554(PR), 0601-0554(PR); where 
(PR) signifies the source that was used as the 
phase reference. The first three pairs contain all of 
the information about the relative source positions 
among the three sources with respect to 0556-1-238, 
and it is these data that were used in order to be 
directly comparable to the VLBA observing strat- 
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egy. However, the observations for the last two 
pairs are independent, and provide an estimate of 
the errors associated with the source positions and 
structures. 

For all VLBA 23 and 43 GHz sessions, the cor- 
related data were sampled every 1 second, with a 
total bandwidth of 64 MHz (2-bit sampling; except 
for 20 and 22 Dec 2008 with 1-bit sampling) con- 
tained in four contiguous frequency bands, each 
with sixteen 1-MHz channels. The center frequen- 
cies were at 22.45 GHz and 43.41 GHz. At 23 GHz 
the amplitude and phase bandpass over the fre- 
quency range were determined from the observa- 
tions of 0556+238, and data that were associated 
with mal-functioning antennas also were flagged. 
At 43 GHz, we inserted two observations of the 
strong ICRF2 source, 0851-1-202, in order to de- 
termine the bandpass amplitude and phase at this 
frequency. For the final VLBA session at 8 GHz of 
one hour, the same frequency and correlator con- 
figuration were used, and the average frequency 
was 8.64 GHz. 

Updates of telescope positions, improved EOF, 
phase compensation for the parallactic angle 
changes over the session, and corrections from 
the monitored antenna electronic gain and phase 
variations were applied to the correlated data. For 
the VERA observations, the phase and delay dif- 
ference between the dual-beams were measured 
using an internally-generated broad-band noise 
source, and these corrections were applied. The 
offsets are large, although relatively stable with 
time, and are critical corrections needed for pre- 
cise phase referencing. 

The estimates of the tropospheric delay for 
the VERA data were obtained from GPS mea- 
surements at each antenna. These are accu- 
rate to about Icni, except at the more hu- 
mi d locations at Ogasawra a nd Ishigaki-jima 



delay res idual over the observation per iod of 40 



minutes ( Mioduszewski fc Kogan I l2009f ). These 



(|Honma. Tamura fc Reid 1120081 ). In order to im 



prove the tropospheric delay model used in the 
VLBA correlator, analysis of the DELZN obser- 
vations, taken near the mid-point of each session, 
determin ed the zenith-pat h delay offset for each 



antenna ( Reid et al 20091 ). The frequency con- 



figuration consisted of four bands each with 16 
1-MHz channels, centered at 23.02 23.57, 23.29 
and 23.46 GHz. The antenna-based group delays 
were measured, and then the data for each antenna 
were fit to a linear clock delay and a zenith-path 



results are not as accurate as the GPS-measured 
tropospheric delays, but decreased the systematic 
zenith-path delay model errors for each VLBA an- 
tennas from about 4 cm, to about 1 cm. At 43 
GHz, zenith-path delay error of 1 cm for a source 
pair two degrees apart give a systematic phase 
error of 18° at the zenith and 30° at a source el- 
evation of 50°. When this error is averaged over 
all antennas, the systematic position shift from 
this uncertainty is estimated to be about 0.01 mas 
(reference to meeting in Socorro). 

These data were then ready for phase-referencing 
calibrations. 

3. Phase Reference Calibrations 

3.1. The Primary CaUbrator 0556+238 

The choice of 0556+238 as the primary refer- 
ence calibrator is discussed in § 12.21 Its ICRF 
history suggested that it is relatively stable over 
time, and was included as a defining sources in 
the ICRF2 catalog since its position varied by less 
than 0.3 mas in both alpha and delta at 8.6 GHz 
over several decades. Although > 80% of the emis- 
sion is contained in a small, unresolved radio com- 
ponent, images of the source were made for each 
frequency and session in order to determine the 
small structure phase. The source peak brightness 
is assumed to be located at the position given in 
Table [T] at both frequencies and for all sessions, 
and this defines the relative position grid for the 
entire experiment. 

It will turn out that the assumption — that the 
location of the peak intensity for this source is 
fixed over the year period of this experiment — is 
incorrect at the level of nearly 0.2 mas! This is sig- 
nificantly less than the image resolution, but the 
source position change with time was easily deter- 
mined from the anomalous motions of the other 
sources. 

3.2. VERA Phase Referencing 

Since both phase calibrator and target are ob- 
served simultaneously with VERA observations, 
phase referencing requires the transferal of the 
measured antenna-based temporal phase and gain 
determined for the calibrator directly to the tar- 
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get. For each 8-min scan on a source pair, the 
AIPS task FRING determined the residual phase, 
and phase derivative with time (rate) and fre- 
quency (del ay) for each on e-minute interval dur- 
ing a scan ()Greisen 2003 ). An example of the 
delay solutions for the a VERA session is shown 
in Fig. [2j The delay values were slowly chang- 
ing with time with an rms scatter over a seven- 
minute scan of about 1.0 usee, caused in roughly 
equal amounts by the limited signal-to-noise over 
one minute and by fast changes in the tropospheric 
refraction. Over the bandwidth of 128 MHz, this 
delay scatter corresponds to a phase change of 30° 
across the band, which causes a decorrelation of 
only about 4%. Delay solutions that were outside 
of the range or with significantly increased scatter 
were flagged. These data usually correspond to 
periods of poor weather conditions. 

After the phase calibrations tables (delay, 
phase-rate and phase) were then applied to the 
phase calibrator. The AIPS task CALIB was 
then run on this data, using the input model of 
05564-238, to determine the gain variations for 
each antenna. Variations of about 10 to 20% over 
the 6.5-hour session were typical. The phase cal- 
ibration tables and the antenna-gain corrections 
found for the phase calibrator were then applied 
to the target source in order to obtain calibrated 
data. Since both sources were observed simultane- 
ously, no temporal interpolation of the corrections 
found for the calibrators were needed to apply to 
the target. 

3.3. VLBA Phase Referencing 

From the "DELZN"40-min observations near 
the mid-point of each session (§ 12. 3p the residual 
zenith-path delay offset and clock delay for each 
antenna were obtained. These corrections were ap- 
plied to all data (23 and 43 GHz) for the session to 
remove the relatively large offset troposphere and 
clock delay error. 

The VLBA phase referencing analysis, using 
the observations of 0556-1-238, were similar to 
those described above for the VERA session. 
These calibrations were done independently for 
23 and 43 GHz. Because of the excellent stabil- 
ity of the VLBA and the relatively dry locations 
of most of the antennas, the AIPS task FRING 
determined only the delay and phase for each 20- 
sec scan of the calibrator. The delays were then 



smoothed over a time scale of one-hour to decrease 
the variations from the signal-to-noise. 

The antenna-based phase variations for the cal- 
ibrator over a session varied considerably, de- 
pending primarily on the weather conditions. 
Examples of these phase variations for several 
VLBA antennas at 43 GHz observations on 08 
June 2009 over a three hour period are shown in 
Fig. [21 In order to maintain the necessary tempo- 
ral phase coherence, the observed antenna-based 
phase must be accurately interpolated between 
adjacent 0556-f238 scans, 55 seconds apart, in 
order to transfer them to the intervening target 
scans. In Fig.[3]periods of poor phase stability are 
obvious (SC 14:30 to 14:50, FD after 16:20, NL 
after 16:50), and the data for the targets during 
these periods were flagged. As with the VERA ob- 
servations, these phase calibrations were applied 
to 0556-1-238 and then the average antenna-based 
gain calibrations were determined using the AIPS 
task CALIB. The delays, phase and gains from 
0556-1-238 are then interpolated between scans in 
order to apply them to the target source. The 
target sources were now ready for imaging. 

4. Radio Source Images 

4.1. Imaging Techniques and Astrometry 
Analysis 

Images for each source, frequency, and array 
were then made from the VERA and VLBA data. 
All image sizes were 1024 x 1024 pixels, with a pixel 
separation of 0.01 mas. The data were weighted 
by the inverse-square of the antenna noise, and an 
additional weighting (called robust=0) of the data 
as a function of was data resolution emphasized 
the fine-scale structure without significant loss of 
signal-to-noise. The dirty images were cleaned us- 
ing small boxes around the areas of known emis- 
sion. Because of the substantial differences in the 
(u-v) coverage among the individual VLBA ses- 
sions and the VERA sessions due to flagged data 
or missing antennas, we used a conservative res- 
olution for the restoring clean beam in order to 
produce images with the same resolution. The 
following resolutions (north/south by east /west 
Full- Width at Half-Power [FWHP] in mas) were 
adopted: for the VLBA at 23 GHz, the resolution 
was 1.0 X 0.5; for the VLBA at 43 GHz, the res- 
olutions was 0.6 X 0.3; for the VLBA at 8.6 GHz 
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the resolution was 3.0 x 1.5; and for the VERA 
23 GHz, the resolution was 1.0 x 1.0. Although 
many radio components in these diagrams appear 
to have an angular size that is similar to that of 
the beam, even a small broadening is indicative of 
internal structure, and both imaging and visibility 
model fitting can determine more precise compo- 
nent parameters and sizes (see § 16.21) . 

The target VLBA images (0547+234, 0554+242, 
0601+245) are affected by phase noise introduced 
by the phase-referencing calibration interpolation 
errors between calibrator scans, that are implied 
in Fig. [3l The VERA observations are not affected 
by interpolation errors. Also, angular decorrela- 
tion (the tropospheric phase refraction between 
calibrator and target are somewhat different be- 
cause of their separation in the sky) corrupts the 
target phases for both the VERA and VLBA ob- 
servations. Unless these phase errors are larger 
than about one radian (all such data are flagged 
before images are made), the phase-referenced 
(PR) images for all sources are on the reference 
frame tied to the apriori position of the phase- 
reference calibrator to an accuracy that can be a 
good as 0.02 ma s over a few degrees on the sky 
(|Fomalont I [20051 ) 

The image quality of the target sources can 
be improved using self-calibration techniques. 
By starting with the PR image as the estimate 
of the source structure and position, the self- 
calibration algorithm then determines improved 
antenna-based phases that are consistent with the 
PR image. The improved antenna-based phases 
are determined every one minute, and corrected 
visibility data are obtained. This technique is 
robust because it utilizes data from all baselines 
(45 maximum for the VLBA) in order to obtain 
antenna-based p hase corrections (9 ma ximum for 
the VLBA). See lCornwell fc Fomalonti (,1999) for 
more information about self-calibration and its 
stability. 

The self-calibrated (SC) images have better 
quality than the original PR images. The peak 
brightness of the stronger components increase 
from a few percent to 20 percent because the cor- 
rected visibility data have smaller phase errors, 
and the extended emission and fainter components 
are somewhat better defined. The improved im- 
age quality aids in determining the angular size 
and frequency spectrum of various components in 



a source. 

Tests from these VLBA and VERA observa- 
tions (as well as from many other VLBI sessions) 
conclusively show that the self-calibration algo- 
rithm, when starting with the phase-referenced 
image, does not change significantly the position 
of the radio emission in the improved images. In 
Fig. in the PR-image and the SC-image are shown 
for 0547+234 at 43 GHz from the 08-Jun-2009 ses- 
sion. This example was chosen because the sep- 
aration between 0556+238 and 0547+234 is the 
largest of the pairs, about two degrees, where 
residual phase errors can become as large as a ra- 
dian at 43 GHz. For this particular example, the 
difference in the location of the peak brightness 
for PR and the SC image is 0.003 mas east/west 
and 0.015 mas north/south. The position differ- 
ence of the secondary peak is 0.02 and 0.07 mas, 
well within the uncertainty of the location of this 
faint and blended component. For other source 
pairs and at 23 GHz frequencies, the differences 
between a PR and an SC image will be less than 
that shown in this example. Hence, in this pa- 
per, the position of all components and the image 
displays will be obtained from the self-calibrated 
data. 

Several methods can be used to determine a 
"position" of a component from a clean image 
including: 1) Measure the location of its peak 
brightness. This position determination is resolu- 
tion dependent, and is subject to the blending of 
close components. This uncertainty can be some- 
what lessened if the same resolution is used for all 
sessions, even though this degrades the resolution 
for some of them. 2) Fit the emission on an im- 
age with a small number of Gaussian-components 
that are needed to describe the emission prop- 
erties. Each component intensity, position and 
angular size, with error estimates, are then ob- 
tained. The angular size measured from the image 
can be deconvolved (since the images resolution is 
known) in order to obtain an estimate of the true 
sky angular size. 3) Fit the calibrated visibility 
data (self-calibrated) directly with a small num- 
ber of Gaussian-components. The initial guess of- 
ten comes from a clean image of the source. Such 
modeling of the visibility data directly, with no re- 
course to the image and its resolution limits, can 
often determine a component's true angular size 
or more complex structure more accurately than 
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by image fitting techniques. 

Method (2) was used for most of the component 
and source analysis. When it is necessary to deter- 
mine the internal structure of some of the compact 
components, method (3) was used in order to ob- 
tain as high a resolution as the data signal-to-noise 
would permit. 

4.2. The Morphology of Each Source 

The VLBA images in the form of contour im- 
ages at 23 GHz and 43 GHz for five sessions are 
given for the phase reference source, 0556-f238, 
and the three target sources, 0547-1-234, 0554-1-242 
and 0601-f245 in Fig. [3 Fig. El Fig.H and Fig.H 
respectively. Although there are six VLBA ob- 
serving sessions, the first contour plot on the left 
shows the average of the images for the 20-Dec- 
2008 and 22-Dec-2008 sessionfl AU images have 
the identical contour levels as a percentage of the 
peak intensity. The quality of the images at 43 
GHz on 12-Dec-2009 and 21-Dec-2008 (the aver- 
age of two sessions) is somewhat degraded because 
of inclement weather at several antennas. This se- 
ries of images show the good consistency of the 
VLBA images at both frequencies, and that none 
of the sources changed substantially over the year 
period. The components for each source have been 
given a number designation. The component '0' is 
thought to be the core, and will be discussed in 
detail in § 15.11 

The images and structure for 0556-1-238 have 
been discussed in § 13.11 in the decision to use this 
source as the main phase reference for the exper- 
iments. It is nearly a point source, although the 
last session at 43 GHz shows a slight extension 
to the south-west. Deeper imaging of this sources 
will show the surprising evolution of this source 
structure over the year observation in §6.21 

The source 0547-1-234 is dominated by a com- 
pact core, displaced about 0.5 mas north-west of 
the phase center, shown by the '+'. A faint com- 



^The intent of these first two sessions, only two days apart, 
were to determine the approximate positional accuracy that 
could be obtained from these observations. A lower sensi- 
tivity set-up was used and the 43 GHz data were not well- 
calibrated. Improvement in the observing strategy were 
then made for the subsequent VLBA sessions. The average 
of these two sessions are comparable in quality with the 
later sessions, so the results are included and are given the 
effective observing date of 21-Dec-2008. 



ponent about 1.0 mas to the west is somewhat 
variable and brightest on the 08-Jun-2009 session. 
The source 0554-1-242 has two components, clearly 
shown in the 43 GHz images. The phase center is 
located near the component to the north-east, but 
the south-west component appears more compact 
and is possibly the radio core. Finally, the source 
0601-1-245 has a complex structure that is resolved 
into three major components at 43 GHz. The pos- 
sible location of the radio core in this source is 
unclear with this level of display. 

4.3. The Frequency Dependence and 
VLB A/ VERA Comparisons 

Most of the analysis in the paper will be de- 
rived from the VLBA observations since these were 
the highest resolution and covered the frequency 
range from 8.6 GHz to 43 GHz. The most ac- 
curate VERA session was on 17-Jan-2009 when 
the weather at the four antennas was excellent. 
For the other VERA sessions one antenna was not 
available for a substantial part of the 6.5-hour ob- 
servation, primarily because of inclement weather. 

The comparison among the VLBA images at 
8.6, 23 and 43 GHz with the VERA images at 
23 GHz are shown in Fig. [9l The importance 
of high resolution is demonstrated in these com- 
parisons. For example, the VERA and VLBA 23 
GHz images show similar structures and position 
for the four sources, but some essential features 
of the source structure are more apparent with 
the better VLBA resolution, especially at 43 GHz. 
The VERA images, although not used heavily in 
the following analysis, were consistent in structure 
and position with the VLBA images, lending more 
confidence to the overall conclusions. The non- 
contemporaneous aspect of the images shown in 
this figure are not a complicating factor in their 
comparison because of the relative slow pace of 
evolution for the sources. 

The VLBA resolution at 8.6 GHz, the same fre- 
quency used for the ICRF, barely hints at the in- 
ner structure of these four sources. Clearly, with- 
out the 43 GHz VLBA images at the highest res- 
olution, interpretation of the source structure and 
evolution would be ambiguous — as will be docu- 
mented in later sections. 
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5. Analysis of the Radio Source Structure 
and Evolution 

5.1. Finding the Radio Core in each Source 

The major goals of the observations of these 
four sources were 1) to determine the location of 
the radio core for each source, 2) to determine 
their position stability, and 3) to assess the af- 
fect of structure on the ICRF2 positions. Finding 
the core components for each of the four sources 
involved somewhat different combination of most 
of the following properties: (1) The most compact 
component with 0.3-mas resolution; (2) The com- 
ponent with the flattest spectral index; (3) The 
component at one end of the emission region; (4) 
The component that is stationary with respect to 
the radio cores in other sources. The use of these 
four criteria will become clear in the following dis- 
cussions. 

The four sources that have been chosen, al- 
though relatively compact, according to the 
ICRF2 catalog, have significantly different struc- 
tural properties with sub-mas resolution imaging. 
The two sources that appear to be core-dominated 
are 0547-1-234 and 0556-1-238, which was used as 
the phase calibrator. On the other hand, the other 
two sources have sufficiently complex structure 
so that the location of the radio core within the 
emission is not obvious from the images in Fig. [7] 
and Fig. |51 The multi-frequency, multi-session 
observations of these two sources were needed to 
"find" their radio core. 

5.2. The Structure and Evolution 
of 0556+238 

The source 0556-1-238 is nearly point-like at all 
frequencies and sessions, which is the main reason 
that it was used as the phase reference calibrator. 
The last session at 43 GHz shows a faint extension 
of the emission to the south-west. The VLBA im- 
ages were fit to a single elliptical Gaussian com- 
ponent, labeled in Fig. [H and its parameters 
are listed in Tabled Column 1 gives the session 
data and the next five columns give the compo- 
nent total flux density, peak flux density and the 
estimated error of the peak, and the offset of the 
component X (east) and Y (north) with error es- 
timates, all at 43 GHz. The next five columns 
give the same intensity and position estimates at 



23 GHz. The ratio of the total flux density to 
the peak flux density, slightly above unity, can be 
caused by the remaining small phase errors (decor- 
relation) or from a non-zero angular size of the 
component which corresponds to a full width at 
half-power diameter of about 0.2 mas. Over the 
experimental period of about one year, the flux 
density of the radio core increased 46% at 43 GHz 
and 18% at 23 GHz, with the spectral index, a 
(where flux density oc i^"), changing from —1.2 to 
-0.8. 

5.3. The Structure and Evolution 
of 0547+234 

The source 0547+234 is dominated by a com- 
pact component toward the east. At 43 GHz, the 
source structure was fit with an elliptical Gaus- 
sian component for component (clearly the radio 
core) , and with a point component for component 
1 (it is too weak for an angular size to be deter- 
mined). At 23 GHz the bright component radio 
core and the fainter component 1 were each fit 
with point source. An additional component 2, 
located about 1.2 mas to the west, was added to 
the source model in order to fit the more extended 
emission that is seen in Fig. [6] The parameters as- 
sociated with components 1 and 2 are relatively 
inaccurate since they are not particularly discrete 
features and appear to change quickly with time, 
but are included in the fit in order to remove any 
interaction of their parameters with that of the 
radio core. These parameters have not been in- 
cluded Table ID Observations of the source at 23 
GHz in Aug 2005 show a similar structure with a 
dominant c ore and a faint compo nent about 1 mas 



to the east ( Chariot et al 2010f) . 



The parameters for the radio core for the five 
VLBA sessions at 43 GHz and 23 GHz are shown 
in Table m with the same format as that for Ta- 
ble [3] The component fiux density reaches a 
minimum for the 08-Jun-2009 session at both 43 
and 23 GHz, when components 1 and 2 are near 
their maximum, about 15% of the flux density of 
the radio core, and the component spectral index 
is about -0.6. It is offset from the phase center 
by about X=-0.18 and Y=+0.21 mas. The most 
surprising property of component is its proper 
motion to the south at 23 GHz and 43 GHz, ap- 
proximately [ly = 0.12 mas/yr.! Slight motion to 
the west may also occur. The cause of this appar- 
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ent motion is discussed below. 

5.4. The Structure and Evolution 
of 0554+242 

The contour plots of 0554+242 in Fig. [7] show 
two major components. The emission was fit with 
a two-component elliptical Gaussian model for all 
sessions and at 43 and 23 GHz, and the parame- 
ters are shown in Table [5] Component is only 
slightly resolved since the total to peak flux den- 
sity is about 1.1, corresponding to an angular size 
< 0.2 mas. There is no significant variability over 
the year experiment. The X-position of this com- 
ponent shows a small motion to the west, but, as 
with 0547-1-234, this component is moving to the 
south by about 0.09 mas over the year. 

Component 1 is clearly extended with the total 
flux density to peak flux density ratio about 1.9 
at 43 GHz and 1.3 at 23 GHz. This ratio corre- 
sponds to an angular size of about 0.8 mas at both 
frequencies. The spectral index of the two compo- 
nents are significantly different. For component 0, 
a = —0.7; for component 1, a — —1.2. It is, thus, 
extremely likely that the radio core of 0554-1-242 is 
component because it is more compact and has 
a flatter spectrum than component 1. 

The relative proper motion between the two 
components (component 1 - component is — 
-1-0.018 ± 0.012 mas/yr and fiy = -f 0.010 ± 0.015 
mas/yr, or /i = 0.020 ± 0.013 mas/yr in position 
angle 56°; but it is only marginally detected. The 
component separation in June 2009 is 0.782 mas in 
position angle 50° which is fortuitously in the same 
direction as the possible motion. Observations at 
23 GHz with the VLBA in 2005 August show a 
separation of the two components of 0.75 mas in 
position angle 54° with flux densities .29 (NE) 
and 0.10 (SW-core) (jCharlot et alll2010[ ). Earlier 
observations at 8.6 GHz from the 1990's also sug- 
gest a similar separation of two major components. 
Thus, the two components are barely if at all sep- 
arating, with an ejection time that was at least 20 
years in the past. 



of the source and an extended component to the 
west (see Fig. [8]). The location of a radio core, per- 
haps corresponding to one of these components, 
cannot be easily deduced from this display with- 
out deeper analysis. 

Table [6] lists the parameters for the three com- 
ponents. The two components to the east were 
modeled as point sources (when two components 
overlap, determination of their individual angular 
sizes is not reliable), but the western component 
was fit to an elliptical Gaussian component, and it 
is clearly extended. The component flux densities 
are relatively stable at 43 GHz, but appear to be 
slowly decreasing at 23 GHz. Component has 
the flattest spectral index with with a = —0.7, 
component 1 has a steep index of a = —1.7, and 
the western, extended component 2 has a — —0.8. 

Component shows a similar proper motion as 
the radio cores in 0547+234 and 0554+242 with 
At^ = -0.029 ±0.015 mas/yr, fiy = -0.064 + 0.025 
mas/yr. Again, this anomalous proper motion is 
discussed in the next sections. This component is 
likely the radio core since it also has the flattest 
spectral index, is compact, and is at one end of 
the radio emission. 

Component 2 is moving with = —0.183 ± 
0.025 mas/yr, fiy = -0.064 ± 0.033 mas/yr, or 
fj. = 0.193 + 0.029 mas/yr in position angle -109°. 
The separation of component 2 from is X = 
-1.231 + 0.022, Y = 0.236 + 0.040, 1.253 mas 
in position angle —101°, which extrapolates to 
an ejection time of component 2 from component 
about 6.5 years ago. Observations at 23 GHz 
with the VLBA in 2005 August shows the pres- 



ence of only one component 



with these measurements (Chariot et al 



which is c onsist ent 
[2010h . 



The motion of component 1 away from is 
= -0.036 + 0.020 mas/yr, fiy = -0.029 + 0.030 
mas/yr, and is formally not detected. Although 
the separation is in the sense that component 1 
is moving away from 0, the blending of the two 
components, even at 43 GHz, makes determining 
their separation uncertain. 



5.5. The Structure and Evolution of 
0601+245 

The 43 GHz images of 0601+245 suggest that 
the source is composed of three major components. 
There are two close components in the eastern part 
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6. Anomalous Position Behavior 
for 0556+238 

6.1. Consistency of Target Core Apparent 
Proper Motions 

The radio cores in 0547+234, 0554+242 and 
0601+245 have been tentatively identified. For 
0547+234 this association was straight-forward 
since the emission is dominated by a single com- 
ponent. For 0554+242 and 0601+245 the 23 GHz 
and 43 GHz high resolution images were both 
needed to determine the core location: it was the 
most compact component, and with the flattest 
spectral index. What is surprising is that cores 
in the three target sources have nearly the same 
proper motion, mainly to the south! 

To investigate these proper motions in more de- 
tail, the positions measured for the core compo- 
nents from Tables Hl5] and [SI are plotted in Fig.lTUl 
Both the 43 GHz and 23 GHz core positions were 
used, although there could be a small separation 
between them, but this should not bias the proper 
motion fit. The relative offset of the cores from 
their phase center has been removed, so the mo- 
tions of the three sources could be overlapped. 
These core offsets were: 0547+234 (X=-0.171, 
Y=+0.194) mas; 0554+242 (X=-0.454, Y=-0.429) 
mas; 0601+245 (X=+0.772, Y=+0.270) mas. The 
first session and the last session points show larger 
scatter than the middle three sessions, as sug- 
gested from the contour plots for these two ses- 
sions in Fig. [SJIHI 

The proper motion fit to the core positions of 
the three target sources are fXx = —0.032 ± 0.013 
mas/yr, ^ly = —0.093+0.017 mas/yr, and aU three 
cores are consistent with this motion. Thus, the 
relative separation among them remained constant 
to the measurement error of about 0.02 mas, as 
indicated by the scatter of the points around the 
best fit proper motion line in Fig. 1101 The similar 
core motion for the three target sources strongly 
implies that the emission peak of 0556+238, which 
defines the location of its apriori position, is mov- 
ing with respect to the other three cores. In the 
next section, super-resolution imaging and mod- 
eling of the emission from 0556+238, will show 
how internal changes in this source structure have 
mimicked a proper motion during the experimen- 
tal period. 



6.2. Super-resolution Imaging of 0556+238 

Although 0556+238 is nearly a point source 
with the nominal resolutions at 23 GHz and 43 
GHz, the source intensity does vary significantly, 
with a change in spectral index, over the year ex- 
periment. Hence, there is some indication of in- 
ternal changes, even within the resolution of the 
43 GHz observations. The source has been well- 
observed for decades in the ICRF-programs with 
little indication of extended structure. The de- 
rived positions from the hundreds of ICRF exper- 
iments are stable at the 0.2 mas level, although 
this is the typical position error for any one ICRF 
session. It has a structure index of 1, meaning 
that it has no significant structure outside of the 
core. It is one of the 295 ICRF2 sources used to 
define the celestial reference frame. 

Errors in the data processing and analysis that 
could produce the apparent motion of 0556+238 
are unlikely. Unexpected astrometric and tropo- 
spheric errors should produce errors that are func- 
tions of the calibrator target separation-and orien- 
tation, rather than a constant offset for all targets. 
An error in the data processing and the correlator 
model in which an anomalous proper motion was 
inserted in the position of 0556+238 over the year 
experiment has not occurred. 

The 43 GHz and 23 GHz data for 0556+238, 
after editing and self-calibration in AIPS, were 
transferred to the Caltech inter ferometric pack- 
age, Difmap (jPerson et al 1 1 19941 ). The visibility 
data were then fit with one or two Gaussian com- 
ponents. For the first two sessions when the source 
was small in angular size, only one extended Gaus- 
sian was fit; but for the three final sessions, when 
the source became more extended, a double point 
model was used. These parameters are shown in 
Table [7] and the images, with a restoring beam of 
0.2 X 0.1 mas at both frequencies, are displayed 
in Fig. [TTJ This resolution is about a factor of 2 
smaller than the conservative one used for previ- 
ous image display at 43 GHz and a factor of 3 at 
23 GHz in Fig. [3 

The NE component is the brighter than the SW 
component by a factor of three at 43 GHz and five 
at 23 GHz, with the spectral index of the NE and 
SW components a ~ —0.8 and a w 0.0, respec- 
tively. At both frequencies the total fiux density of 
the source at both frequencies increased by about 
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50% from Dec 2008 to Dec 2009. The '+' in the 
contour images of 0556+238 in Fig. [TT] shows the 
location of the assumed apriori position for the 
source. The location is determined by the self- 
calibration algorithm and is usually within 0.005 
mas of the peak brightness. Since the NE compo- 
nent is much brighter than the SW component, the 
apriori phase center is virtually coincident with 
the NE component. 

The source contains two components that are 
moving apart. From the model fits the component 
separations are 0.15, 0.19, 0.24 mas for the June, 
Sep and Dec 2009 sessions, with good agreement 
between the 23 and the 43 GHz results. Using the 
ratio of peak to integrated flux densities for the 
first two sessions as a measure of the separation of 
the two components, we obtain about 0.10 mas for 
both sessions. Thus, to the accuracy of these fits, 
the proper motion between the two components is 
about 0.14 ± 0.02 mas/yr in position angle 150 ± 
15°. For comparison, the average proper motion of 
the three cores in the target sources (See Fig. [TO)) 
is 0.10±0.02 mas/yr in position angle 161 ±12°. If 
these three cores are assumed to be stationary in 
the sky, then a stationary point in the images for 
0556-1-238 over the sessions is shown by the 'box'. 
It is clear that the SW component is not moving 
significantly with respect to the other three target 
radio cores. 

The conclusion is that the south-west compo- 
nent is the stationary radio core of 0556-1-238. The 
brighter component is a new jet component, mov- 
ing to the north-north east with a proper motion 
of about 0.12 mas/yr. It appears fortuitous that 
this motion is quasi-linear because the first ses- 
sion occurred somewhat after the jet ejection, and 
the jet component velocity has remained relatively 
constant over the year period, a typical behavior 
of many jet compo nents just after their formation 



(|Lister et alll2009l ). 



6.3. Super-resolution Imaging of Other 
Sources 

The visibility data were also analyzed for the 
other sources, using similar super-resolution tech- 
niques that were described above. The cores for 
the three target sources were compact for all of 
the sessions, and were sufficiently far from the 
additional extended emission, so that their posi- 
tion could be unambiguously determined to an ac- 



curacy of 0.02 mas. The unusual circumstances 
affecting the interpretation of the structure of 
0556-1-238 — the ejection of a new component from 
the core during the experiment period, and the 
intensity of this component significantly greater 
than that of the core — did not occur for the other 
three sources. 

7. Comparison with the ICRF2 and the 
8.6 GHz Positions 

The Jan 2010 VLBA session observed the four 
sources in phase-reference mode at 8.6 GHz for 
two main reasons: First, to obtain the 8.6 GHz 
structures and positions that are relatively co- 
temporaneous to the 23 and 43 GHz images in 
Dec 2009, so that time variability is not a serious 
problem for their comparison. Second, to assess 
the ICRF2 positions with the detailed structure 
information obtained for the four sources. 

7.1. ICRF2 Registration of Images using 
0556+238 

An overlay of 0556+238 of the 43-GHz image 
in Dec 2009 and the 8.6-GHz in Jan 2010 image is 
shown in Fig. [12] The observations were taken 45 
days apart, so the slight variability or motion in 
the source should not be a contributing factor in 
the image differences. 

However, the 43 GHz and 8.6 GHz grid frames 
could be displaced. This is unlikely because the 
assumed phase center (same for all three frequen- 
cies) is virtually located at the position of the 
peak brightness, and the NE component is much 
brighter than the SW component at 8.6, 23 and 
43 GHz. Thus, the Dec 2009 23 and 43 GHz 
images and Jan 2010 8.6 GHz images are on the 
same reference grid to better than 0.1 mas. This 
means that the comparison of the 8.6 GHz and 43 
GHz images for the other sources will also be ac- 
curately registered since these are all tied to that 
of 0556+238, the phase reference calibrator. 

7.2. Where is the ICRF2 Position of 
0556+238? 

The ICRF2 position of 0556+238, one of 
the 295 sources that are used to define the 
ICRF quasi-inertial reference frame, is a = 
05''59™32.03313165",^ = 23°53'53.9267683" with 
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an estimated error of 0.04 mas in both coordi- 
nates. This improved position is 0.026 mas east 
and 0.068 mas north of the assumed phase center 
(see Table that was taken from the previous 
ICRF catalog, and the ICRF2 position is 0.07 mas 
displaced in the same direction as that between 
the NE and SW components! 

An estimate of the evolution of the structure 
of this source structure can be made from VLBA 
observations at 23 and 43 GHz in Dec 2002, 
Sep 2003 and from 23 GHz and 8.6 GHz from 
Aug 2005 (|Lanvi et al ll2010HCharlot et al Il2010[ ). 
Over-resolved images of the source for these ob- 
servations show two components that were sepa- 
rated between 0.1 and 0.3 mas in a position angle 
roughly —135°. The relative flux density between 
the two components was not fixed, but the south- 
ern component had a flatter spectral index in two 
of the three experiments, which supports the con- 
tention that the radio core is indeed the south-west 
component. The time-scale of the formation and 
motion this jet component cannot be ascertained. 

The location of the ICRF2 position of 0556+238 
could lie anywhere between the NE and SW com- 
ponents, although at 8.6 GHz, the NE component 
is probably the stronger of the two components 
and the ICRF2 position may favor this position. 
The error cross in Fig. [12] is a conservative esti- 
mate of the range of the location of the ICRF2 
position in Dec 2009 to Jan 2010. This error box 
will be transferred to the target sources in order 
to compare their 23 and 43 GHz images with the 
8.6 GHz images and with the ICRF2 position. 

7.3. 0547+234 at 8.6 GHz 

A comparison of the 23 GHz Jun 2009 and and 
8.6 GHz Jan 2010 VLBA images is given in Fig.HSl 
The June 2010 session at 23 GHz was used since 
it shows best the extended emission to the west. 
(The motion of the radio core between June and 
Dec 2009 was less than 0.1 mas which does not 
affect the registration of the two images.) Two 
anomalies are apparent. First the ICRF2 position 
is significantly south of the peak emission at all 
frequencies. This source has been observed in only 
about 23 ICRF sessions over the years, and the 
formal ICRF2 position error is 0.3 mas. Hence, the 
offset of the ICRF2 position could be consistent 
with the images. 



The second anomaly, the displacement of the 
peak brightness at 8.6 GHz and 43 GHz (and all of 
the high-frequency images), is significant. As seen 
in Fig. [SI the faint emission to the west comes 
and goes over 2009, at both 23 and 43 GHz. A 
radio image at 8.6 GHz in Aug 2005 shows even 
more extended structure to the west of the peak 
(Ihttp: / / rorf.usno.navy.mil/RRFID / ) than that in 
Jan 2010. An explanation for this discrepancy is 
that there is significant diffuse and variable radio 
emission emanating north and west of the radio 
core. This emission may be optically thick within 
1 mas of the core and, is thus, barely detected 
at 23 and 43 GHz. Thus, the change of position 
of the core between 43 GHz and 23 GHz may be 
the so-called frequency depe ndent core-shifts, se e 
for some other radio sources ( Kovalev et al |[2009l) . 
which is typically about 0.2 mas between 23 GHz 
and 8.6 GHz. This is the only one of the four 
sources which has extended emission near the ra- 
dio core. Nearly all of the emission in the other 
three sources are confined to discrete components 
and such core-shifts are not apparent. 

7.4. 0554+242 at 8.6 GHz 

A comparison of the 43 GHz Dec 2009 and and 
8.6 GHz Jan 2010 VLBA images is given in Fig. [Ml 
The 8.6 GHz image, which shows a slight asym- 
metry to the west, can be modeled with two Gaus- 
sian components that nearly coincide with the 43 
GHz component positions. The NE component at 
8.6 GHz contains 0.73 Jy, compared with 0.10 and 
0.24 Jy at 43 and 23 GHz, respectively. The SW 
component at 8.6 GHz contains 0.13 Jy, compared 
with 0.07 and 0.13 Jy at 43 and 23 GHz, respec- 
tively. Hence the steep spectrum of the NE com- 
ponent and the flat spectrum of the SW compo- 
nent are supported by the 8.6 GHz model. Thus, 
the position (peak brightness) of the radio source 
in Jan 2010 is the NE component which is now 0.8 
mas north east of the core. 

The ICRF2 position of 0554+242 is located be- 
tween the two components, although somewhat 
closer to the NE component, but is still 0.6 mas 
from the core position. Although no significant 
motion between the two components has been de- 
tected over the last ten years, the present ICRF2 
position (generated from more than 20 years of 
data) is consistent with earlier periods when the 
NE component may have been closer to the core. 
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or reflects on the different intensity ratios between 
the two components. 

7.5. 0601+245 at 8.6 GHz 

A comparison of the 43 GHz Dec 2009 and and 
8.6 GHz Jan 2010 VLB A images is given in Fig. [HI 
The source structure is complicated, but analysis 
of the 23 and 43 GHz images suggests that the ra- 
dio core is associated with the component at the 
eastern part of the emission. The peak of the 8.6 
GHz image (hence its position) is near the middle 
component which has a relatively steep spectral in- 
dex and would dominate the emission at the lower 
frequency. The western component can be see in 
the 8.6 GHz image. 

The ICRF2 position is displaced significantly to 
the east, more than 1.0 mas from the radio core, 
and about 0.6 mas east of the present 8.6 GHz po- 
sition. The western component that has a proper 
motion of about 0.19 mas/yr which extrapolates 
to a lifetime of about 6 years if the motion has 
remained the same. There are few good quality 
observations of this source over the last five years 
at 8.6 GHz; however, if this western jet compo- 
nent was stronger in the past, the ICRF2 position 
would then have been dominated by this moving 
component. 

8. Discussion 

The four ICRF2 sources chosen in this study ap- 
peared relatively simple from the 8.6 GHz ICRF 
monitoring over several decades. Their emission 
was dominated by a compact component and the 
position error estimates ranged from 0.05 to 0.15 
mas. The observations at 23 and 43 GHz, reported 
here, show that the situation is much more compli- 
cated, even for these high-quality ICRF2 sources. 

Although this high resolution, accurate astro- 
metric, multi-frequency, multi-epoch set of obser- 
vations covered only four sources, some general 
conclusions about the accuracy of the 8.6 GHz 
ICRF and the properties of radio cores can be in- 
ferred. 

8.1. An Evaluation of the ICRF 8.6 GHz 
Positions and Structures 

A comparison of the ICRF 8.6 GHz positions 
and images with those obtained with the present 



observations for each of the sources can be sum- 
marized as follows: 

0556+238: The structure and evolution of this 
source, one of the 295 defining objects of 
the ICRF2, is not simple even though over 
90% of the emission is contained with an an- 
gular size of 0.5 mas. Many years of ob- 
servations show virtually no extended emis- 
sion at 8.6 GHz and perhaps a slight ex- 
tension at 2.3 GHz to the north and west 



( http: //rorf. usno.iiavy.mil/RRFID /\ . Dur- 
ing the course of this experiment, however, 
a strong jet component moved from 0.1 to 
0.2 mas from the weaker, but optically thick, 
radio core. Since the source was used as the 
phase reference calibrator, it produced ini- 
tially what were anomalous properties of the 
three target sources: their radio cores were 
moving in the sky!. This type of structure 
variation within a region of 0.3 mas for this 
source was corroborated by analyzing sev- 
eral other experiments over the last seven 
years. A lesson learned from the results of 
this source, is that structure changes of an- 
gular scale 20% of the resolution can occur, 
and are particularly injurious in determining 
accurate radio positions. 

0547+234: This source contains about 80% of 
its emission in a compact component be- 
tween 8.6 and 43 GHz. However, the posi- 
tion at 8.6 GHz lies about 0.3 mas of the 
radio core at 23 and 43 GHz. The faint 
emission, observed at all frequencies toward 
the west of the core, may be a diffuse com- 
ponent that becomes progressively optically 
thick at lower frequencies. Unlike the other 
three sources, this extended emission is not 
confined to well-defined components, but is 
more ephemeral in time and space. This ra- 
dio core separation between the 8.6 GHz and 
43 GHz frequencies is similar to that see n 

[iool. 



for other radio cores (Kovalev et al 



However, the other three radio sources con- 
tain most of their emission in well-confined 
components and similar core shifts associ- 
ated with an optically thick diffuse emission 
extending from the core do not occur. 

0554+242: The source structure is relatively 
simple from the 23 GHz and 43 GHz obser- 
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vations, containing a compact, radio core to 
the south-west and a (possibly) slow mov- 
ing jet component to the north-east. The 
present 8.6 GHz location coincident with the 
jet component since it is the brighter compo- 
nent at 8.6 GHz. Thus, the present position 
is displaced currently (Jan 2010) about 0.8 
mas from the core. The cataloged ICRF2 
position (based on 10 to 15 year average po- 
sition ) is located somewhat closer to the ra- 
dio core for two possible reasons: The NE 
components was closer to the core 15 years 
ago, and/or the flux density of the radio core 
was more dominant. 

0601+245: The source is relatively large and 
complex, and even the 8.6 GHz-images show 
an extension to the west. Only with the 
43 GHz VLBA observations could the radio 
core be identified as the most eastern com- 
ponent of the three. The present 8.6 GHz 
position is dominated by the bright, steep- 
spectrum component that is about 0.5 mas 
west of the core. The IGRF2 positions is fur- 
ther displaced by 0.7 mas to the west. Unfor- 
tunately, there are no previous observations 
for which an accurate structure at 8.6 GHz 
be obtained. It is suggested that the western 
component, which is moving 0.2 mas to the 
west, was much brighter ten years ago and 
this component influenced the ICRF2 posi- 
tion much more than presently. 

Several conclusions from the study of the four 
sources are: First, the position of the source at 8.6 
GHz (whether from the VLBA Jan 2010 observa- 
tion or the IGRF2) for these relatively compact 
sources is displaced 0.2 to 1.0 mas from the radio 
core. The displacement is usually in the direction 
of the more extended structure. Thus, the 8.6 GHz 
radio positions, often with an estimated error of 
< 0.1 mas, will often not include the position of 
the radio core. The offset may be caused by the 
optical depth affect of a relatively smooth jet con- 
nected with the radio core (0547-1-234 may be an 
example) and is usually designated as a core-shift. 
However, for the other three sources, the 8.6 GHz 
position is influenced by the strong emission from 
a discrete component that is not the radio core. 

Second, the components that are associated 
with the 8.6 GHz positions may move away from 



the core, and hence the 8.6 GHz positions also 
change with time. This change was noted for 
two of the sources where the 8.6 GHz positions 
measured in Jan 2010 significantly differ from 
those of the ICRF2 that are average positions 
for a source over the last 10 to 15 years. Sim- 
ilar changes (small linear motions) in the posi- 
tion of ICRF sources have been observed for many 
ICRF2 sources over the decades of monitoring 
( Feissel- Vernier. M. 20031) . and the present ob- 



servations suggest that such position changes are 
structure-induced, rather than intrinsic to the ra- 
dio source (i.e. radio core) position. 

8.2. Identifying the Radio Cores 

From analysis of 43 and 23 GHz VLBA phase- 
referenced images of four source observed over a 
year, the location of their radio cores have been as- 
certained. The core component for each of the four 
sources had some or all of the following properties: 
the following properties: (1) The most compact 
component with 0.3-mas resolution; (2) The com- 
ponent with the flattest spectral index; (3) The 
component at one end of the emission region; (4) 
The component that is stationary with respect to 
the radio cores in other sources. 

These criteria could be evaluated by observing 
the sources in the following way: 

43 GHz Resolution of 0.3 mas: The long 
VLBA baselines at 43 GHz were crucial in 
resolving and determining the radio cores. 
The resolution of 0.3 mas is needed to sepa- 
rate close components or to resolve out the 
inner jet. For two of the sources the VLBA 
resolution at 23 GHz was not sufficient. A 
second reason for observing at 43 GHz is 
that the emission, even from the core re- 
gion of most AGN's, is becoming opticall y 
thin (jLanvi et al IbOloHCharlot et al II2OI0I) . 
Thus, at 43 GHz little emission is expected 
to be contained in opaque components that 
produce structure changes with frequency. 
This is supported by the spectral index of 
-0.7 for several of the radio cores, although 
a few had fiat spectral indices. 

23 GHz Observations: A second observing fre- 
quency is needed to determine the spectral 
index of each component, as well as produc- 
ing more accurate images than at 43 GHz 
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(better SNR and less tropospheric contami- 
nation) , albeit with less resolution. The two 
frequencies should not be separated by more 
than a factor of two in frequency in order to 
lessen resolution effects in the determination 
of the component spectral index and posi- 
tion. The spectrum of a component is one 
of the indicators of radio core or non-radio 
core emission. The cores have a spectral in- 
dex > — 0.7 while the other components tend 
to have steeper spectra. 

Multiple Session: Multiple sessions, separated 
by months to a few years, are needed to de- 
termine which components have significant 
motion. This clearly aids in the interpreta- 
tion of the source structure and can be used 
among several sources in order to determine 
which components in several sources are mu- 
tually stationary. 

Phase Referencing: The most accurate method 
to determine the relative positions between 
different sources to about 0.02 mas (e/w) 
and 0.03 mas (n/s) is with the phase refer- 
encing. The gain in positional accuracy is 
because the large tropospheric delay vari- 
ations cancel among closely-space sources. 
Unfortunately, phase referencing cannot be 
used for sources more than several degrees 
apart because of the tropospheric refraction 
that varies linearly with source-pair separa- 
tion. Group-delay referencing can be done 
over a much larger region of the sky, but the 
positional accuracy over a session is likely to 
be about 0.1 to 0.2 mas, limited by the large- 
scale tropospheric refraction. As described 
above, the combination of high resolution 
and the ability to measure the relative po- 
sition of components that were in different 
sources several degrees apart was crucial in 
finding which of the components in the four 
sources had separations that were stable to 
the 0.02 mas level. This was strong evidence 
that they were stationary in the sky since 
there is little reason for the components in 
four independent sources to be precisely co- 
moving. 



8.3. Radio Core Stability 

The measured radio core position stability for 
the four sources is illustrated in Fig. [1^1 The 
plot contains the data points from Fig. [TUl af- 
ter removing the proper motion of the reference 
grid caused by the motion of the strong jet com- 
ponent in 0556+238. The data for the SW- 
component of 0556+238 from Table[71 when blend- 
ing with the jet component was minimal, was 
also included. The rms scatter among the ses- 
sions and the sources is dominated by the tropo- 
spheric refraction residuals of about 4 psec (0.02 
mas). The scatter for the first session (20Dec2008 
+ 22Dec2008 sessions combined) and the last ses- 
sion is some larger than the middle three sessions 
because of lower sensitivity observations in Dec 
2008 and unusually poor weather at many VLBA 
locations in Dec 2009. The normalized is 0.7 
and 1.6 for the n/s and e/s offsets, suggesting that 
the n/s estimated errors are a factor 0.8 too high 
and the e/w estimated errors are a factor 1.3 too 
low. 

There is no significant relative motion for any of 
the radio cores, and the rms scatter per point (one 
radio core for one session) is about 0.02 mas e/w 
and 0.03 mas n/s. These sources are separated by 
an average of 1.5°, so a position stability of 0.02 
mas of about one part in 2.7 x 10* was achieved. 

Radio cores are believed to be located at the 
base of the jet which may be about 2 au from the 
massive object and accretion disk near the nucleus 
of the galaxy. This separation is about 0.003 mas 
for the four sources, assuming a redshift of 1.0. 
Whether the radio core position (and the jet base 
position) are functions of time is not known, and 
one of the goals of these observations are to de- 
termine such changes over time. At present, the 
radio cores appear stable to about 0.02 mas over 
one year. Clearly, further observations of these 
four sources are needed to determine the longer 
term stability of radio cores. Phase referencing of 
other sets of sources are also needed to enlarge the 
sample of ICRF sources in which the location of 
the radio core is known with some assurance. 

8.4. The Source Structure Effect On the 
ICRF 

The ICRF has determined the position of sev- 
eral hundred sources from the location of the emis- 
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sion peak associated with compact radio sources at 
8.6 GHz. The comparison of detailed observations 
of four sources at higher resolution and frequency 
than that of the ICRF showed that the positions 
are often displaced up to 0.5 mas from the radio 
core, and could vary in position by 0.3 mas/yr. 
The ubiquity of these changes, determined from 
only four sources, cannot be securely extrapolated 
to the defining sources in the ICRF2 catalog as 
a whole, but these four sources were among the 
most compact and brightest of the ICRF source 
list. However, the ICRF2 grid averages the foibles 
of each of the 295 sources, so that the average off- 
set or motion to the inertial grid is decreased by 
about 15 from any average source anomaly. This 
suggests a stability of 0.02 mas/yr in the orien- 
tation of the grid, and a possible displacement 
of 0.03 mas in the absolute position of the grid 
should be obtained even if the position of many 
of the sources are displaced from the radio core in 
a similar manner to the four sources studied here. 
Hence, the cumulative effect of source structure in- 
duced changes in the ICRF2 grid are comparable 
to the errors suggested from the residual tropo- 
spheric refraction residuals that affect the current 
best modeling of the radio source positions and 
other related astrometric/geodetic parameters. 

When dealing with individual sources, the es- 
timated errors associated with the ICRF2 posi- 
tions are determined from the scatter of many 24- 
hour session position determinations over several 
decades. Some of the scatter are caused by VLBI 
network problems and troposphere modeling er- 
rors and tend to be random over periods of months 
to years. On the other hand, the apparent position 
offset of the ICRF2 position and the radio core, as 
suggested in this paper, are generally biased in the 
direction toward the inner jet which is remains at 
the same orientation with the core over many years 
for most radio sources. Thus, the estimated ICRF 
position errors reflect the more random changes 
and can be as small as 0.04 mas for sources that 
have many observations. Thus, the offset between 
the apparent 8.6 GHz position and the radio core 
has a systemic component so the offset can be sig- 
nificantly larger than the estimate position error. 



the regions of high optical depth at 8.6 GHz which 
hide the location of the true core, and sufficiently 
high resolution to resolve and recognize the ra- 
dio core from other components. Hence, imaging 
the sources in the ICRF2 defining list at 43 GHz 
(even without phase referencing) is recommended 
in order to determine the internal structure and 
changes that are often hidden with the 8.6 GHz 
observations because of limited resolution and op- 
tical depth effects. These observations would also 
aid in weeded out some of the 295 defining sources 
which have the largest position instability. 

8.5. An ICRF at 23 or 43 GHz 

Phase referencing techniques can determine the 
accurate structure and relative position of sources, 
but only over a small angular region of several de- 
grees. Hence, its extension to the entire sky would 
take a long period of time and the coupling of 
phase referencing results over a few degrees into 
a whole sky reference frame would probably add 
zonal systematic errors that would be hard to re- 
move. 

Recently, observations to determine an all-sky 
catalog at 23 and 43 GHz were initiated. Part of 
the motivation was to establish a catalog of source 
positions for future spacecraft tracking (at 33 
GHz), and part was to determine if improved accu- 
racy could be obtained at 33 GHz co mpared with 
8.6 GHz. The results ar e given in (jLanvi et al 



( Feissel- Vernier. M. 112003 ) 



In order to resolve the radio core and to under- 
stand the evolution of an individual radio source, 
higher frequency observations are needed to avoid 



2OIOI : ICharlot et aTlbOlOl) , with most of the ob- 
servations were made at 23 GHz. The basic ob- 
servable was the group delay, and similar reduc- 
tion and analysis techniques, as that for the ICRF, 
were used. 

Although images were made for each source 
for the twelve observing sessions, they were not 
analyzed as completely as those in this paper. 
Many sources were dominated by a strong com- 
pact component with the source position given as 
the peak location of this component. However, 
some of these components could be resolved into 
sub-components, one of which might be the ra- 
dio core that could be displaced more than 0.1 
mas from the cataloged position. Nevertheless, 
these results strongly suggest that observations at 
23 and/or 43 GHz will likely produce images for 
which a radio core can be identified. The analysis 
of the group delay can then be modified to remove 
the structure effects and determine the position 
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for the radio core component and not simply the 
position of the maximum emission. However, until 
the tropospheric modeling can be significantly im- 
proved, it will be the main limitation to positions 
derived at 32 GHz for any single session, and with 
somewhat minimal gains even with incorporating 
the imaging information. 

8.6. VLBI2010 

The project VLBI2010 (|NieU et al II2OIOI ) is the 
next generation radio array that will be used 
primarily to determine the celestial and terres- 
trial reference frames significantly more accurately 
than the present frame. Two improvements in 
observing technique will be used: First, the tele- 
scopes in the array will be small so that they can 
slew over large areas of sky very quickly. This 
will permit a more accurate determination of the 
tropospheric refraction model over the sky more 
quickly and with more angular dependence. Since 
the troposphere is the major limitation to the 
present ICRF, this alone will produce a significant 
gain in accuracy, perhaps by a factor of two. 

The second improvement will be to observe in 
a continuous band from 2.5 to 15 GHz and with 
sufficient SNR in order to convert the measured 
group delays into unambiguous phase delays. Be- 
cause of the significant ionospheric refraction com- 
ponent and the changing structure and resolution 
of the sources with frequency, the fitting of the 
phase delays from these multi-frequency observa- 
tions will be complicated. In other words, the ba- 
sic observable, the delay (linear phase slope with 
frequency) will be contaminated by the ionosphere 
delay (phase versus inverse frequency) and the 
frequency-dependent source structure phase. This 
phase change will consist of the apparent shift- 
ing of the source position with frequency as well 
as changes of its internal structure. Because the 
highest planned observing frequency is 15 GHz, it 
is likely that there will not sufficient resolution to 
resolve the radio core within the observed struc- 
ture. The multi-frequency images may help to de- 
termine which part of the source has the flattest 
spectral index which could possible identify the 
radio core. 

It is recommended, as it is with the current 8.6 
GHz ICRF, that occasional observations with the 
VLBA be made at 43 GHz in order to assess the 
radio structure at the level of 0.2 mas to better 



determine the true position of the radio core. An- 
other choice may be the VERA-I-KVN array, with 
seven 20-m telescopes, that will begin operation at 
43 GHz by 2012. The use of the EVN is somewhat 
limited in its 43-GHz capabilities at the present 

8.7. The Future of the ICRF 



The ICRF as defined by the ICRF-Ex t2 (iFev et al 



2004), and improved with the ICRF2 (|Fev et al 



2010I ). is the realization of the International Refer- 
ence System using extragalactic sources (quasars) 
to establish a quasi-inertial grid at radio wave- 
lengths. This is the primary frame to which the 
positions of all celestial objects arc referenced. 
The 295 objects which define the ICRF2 sys- 
tem are somewhat sparsely distributed over the 
sky. Since these objects are typically fainter than 
16-mag, the primary optical catalog, Hipparcos, 
is referenced to the radio frame through a few 
bright quasars and secondarily with fainter radio 
stars that are a ccurately tied to the ICRF-grid 
( Lestrade et al l[l998. ). The gaia mission will mea- 
sure the positions of starts as faint at 20-mag 
visual magnitude and will detect over 100,000 
quasars in order to p roduce a large qua sar ref- 
erence frame (LQRF) ([Andrei et alll2009[) . 

Gaia is a scanning instrument that will use 
quasars to determine the parallax zero point, the 
frame spin, and the ICRF frame tie. It is esti- 
mated that the mission accuracy will be better 
than 0.025 mas for stars brighter than 16-mag in 
the visual, but decreasing to 0.2 mas for 20-mag 
objects. With the large number of quasars. Gala's 
measurement should define the inertial reference 
frame axes to 0.0001 to 0.001 mas. Since the ra- 
dio core positions in the ICRF catalog may be in 
error by 1 mas, the frame defined by Gaia may be 
superior to that of the radio. 

The number of precise radio positions has in- 
creased dramatically with the use of the Very 
Large Array (VLA0 and the VLBA (jPetrov et al ' 
20081 ). When combined with the ICRF catalog 



there are over three thousand reference points in 
the sky. These catalogs gives a good densification 
in the northern hem isphere, and is now increas- 
ing m the south (Petrov et alll2009l) . With this 
radio densification, comparison of the Gaia iner- 
tial frame with the radio frame may detect zonal 



•^http: / / www.vla.nrao.edu / astro / calib / manual / 
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errors in either or both frames. 

9. Conclusions 

Using VLBA phase referencing on four close 
ICRF2 sources at 23 and 43 GHz for five ses- 
sions over one year, we have identified the radio 
cores and determined the radio structure of each 
source with 0.3-mas resolution, found the spectral 
index of the radio components, and measured the 
motion of the radio components in each source. 
The accuracy of the relative positions of the four 
sources (through their radio cores) was about 0.02 
mas. The fours cores showed no relative motion 
over the year within the astrometric accuracy. The 
sources contained several other components, some 
of which were moving up to 1.2 mas/yr from the 
radio core. 

With one VLBA 8.6 GHz session and the 
ICRF2 positions for these radio sources, the 
ICRF2 position determination and stability were 
assessed. We found that the ICRF2 position can 
be dominated by a jet component, rather than the 
radio cores that were detected at 43 GHz. Thus, 
the ICRF2 position may be offset from the radio 
core by up to 0.5 mas. Since many of the jet com- 
ponents are moving, the ICRF2 positions will also 
change, perhaps by 0.2 mas/yr, until the jet com- 
ponent becomes less intense than the radio core 
or becomes sufhciently separated (typically by 1.0 
mas). Such motions have been seen in the ICRF 
data base, and all can be attributed to source 
structure. 

Until significantly better tropospheric mod- 
eling of ICRF-type observations can be made, 
it will be the dominant source of error, even 
with these astrometric errors noted in this pa- 
per. Recent VLBA 23 GHz ICRF-type obser- 
vations have obtained similar accuracies to the 
8.6 GHz ICRF2, but more accurate imaging and 
analysis are needed to find the radio cores within 
the brightest components. The VLBI2010 project 
shows promise in increasing the accuracy of the 
reference frame and other astrometric parame- 
ters, but should consider extending their highest 
frequency of 15 GHz to 23 or 43 GHz in order 
to determine the stable radio core location with 
more confidence. 
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Fig. 1. — The location of the four ICRF sources 
used for these observations 
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Fig. 2.— The Delay Solutions for the VERA 
17-June-2009 session. For each 8-min scan of 

0547+234, a delay solution every 1 min was ob- 
tained. Mizusawa is in northern Honshu, Iriki is 
the reference antenna and is located 1000 km wast 
of Tokyo. Ogaswara is located 1000 km south of 
Tokyo, and Ishigaki is located 200 km east of Tai- 
wan. 



Fig. 3.— The Phase Stability at 43 GHz for ob- 
servations on 14-Sep-2009 for the 0556+256 scans 
between UT 14:20 to 17:20 for selected antennas. 
Each point is the average antenna-based phase for 
an 18-scc scan of 0556+238, with LA as the ref- 
erence antenna. Periods of poor phase stability 
between adjacent scans are: nearly all data after 
UT 16:45, SC between 14:30 and 15:00. The large 
phase slope for MK can be followed between cali- 
brator scans. 
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Fig. 4. — A comparison between the phase- 
reference (PR) image (left) and the self-calibrated 
(SC) image (right) for 0547 at 43 GHz on 08-Jun- 
2009. The SC image was obtained by starting with 
the PR image and then determining the residual 
antenna-based phase errors most consistent with 
the PR image. Both images have been cleaned 
and the restoring beam for both is 0.6 x 0.3 mas. 
The contour levels are at -3,3,6,12,24,50,75,99.5% 
of the peak brightness of 0.087 Jy and 0.115 Jy 
for the PR and SC images, respectively. The grid 
lines are separate by 0.5 mas. 
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Fig. 5. The VLB A images at 23 GHz and 43 
GHz for the phase reference source 0556+238: The 
top row shows the 23 GHz contour plots from left 
to right for the sessions: 2008-Dec-21, 2009-Mar- 
16, 2009-Jun-08, 2009-Sep-14 and 2009-Dec-09; re- 
spectively. The bottom row shows the 43 GHz 
contour plots. The resolution at 23 and 43 GHz 
is 1.0 X 0.5 mas and 0.6 x 0.3 mas, respectively. 
The field of view for all images is 3.4 x 2.6 mas, 
and the tick marks are separated by 0.5 mas. The 
'+' shows the location of the phase center for the 
source. All contour levels are at -4, 4, 8, 16, 32, 
50 (bold), 80, 99.5 percent of the peak flux density. 
This component is designated as 0. 




Fig. 6. The VLBA images at 23 GHz and 43 
GHz for the phase reference source 0547+234: The 
top row shows the 23 GHz contour plots from left 
to right for the sessions: 2008-Dec-21, 2009-Mar- 
16, 2009-Jun-08, 2009-Sep-14 and 2009-Dec-09; re- 
spectively. The bottom row shows the 43 GHz 
contour plots. The resolution at 23 and 43 GHz 
is 1.0 X 0.5 mas and 0.6 x 0.3 mas, respectively. 
The field of view for all images is 3.4 x 2.6 mas, 
and the tick marks are separated by 0.5 mas. The 
'+' shows the location of the phase center for the 
source. All contour levels are at -4, 4, 8, 16, 32, 
50(bold), 80, 99.5 percent of the peak flux density. 
The three components contained in the source will 
be designated as 0, 1, 2. 
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Fig. 7. The VLBA images at 23 GHz and 43 
GHz for the phase reference source 0554+242: The 
top row shows the 23 GHz contour plots from left 
to right for the sessions: 2008-Dec-21, 2009-Mar- 
16, 2009-Jun-08, 2009-Sep-14 and 2009-Dec-09; re- 
spectively. The bottom row shows the 43 GHz 
contour plots. The resolution at 23 and 43 GHz 
is 1.0 X 0.5 mas and 0.6 x 0.3 mas, respectively. 
The field of view for all images is 3.4 x 2.6 mas, 
and the tick marks are separated by 0.5 mas. The 
'+' shows the location of the phase center for the 
source. All contour levels are at -4, 4, 8, 16, 32, 
50 (bold), 80, 99.5 percent of the peak flux density. 
The two components contained in the source will 
be designated as 0, 1. 




Fig. 8. The VLBA images at 23 GHz and 43 
GHz for the phase reference source 0601+245: The 
top row shows the 23 GHz contour plots from left 
to right for the sessions: 2008-Dec-21, 2009-Mar- 
16, 2009-Jun-08, 2009-Sep-14 and 2009-Dec-09; re- 
spectively. The bottom row shows the 43 GHz 
contour plots. The resolution at 23 and 43 GHz 
is 1.0 X 0.5 mas and 0.6 x 0.3 mas, respectively. 
The field of view for all images is 3.4 x 2.6 mas, 
and the tick marks are separated by 0.5 mas. The 
'+' shows the location of the phase center for the 
source. All contour levels are at -4, 4, 8, 16, 32, 
50 (bold), 80, 99.5 percent of the peak flux density. 
The three components contained in the source will 
be designated as 0, 1, 2. 
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0556+238 



0547+234 



0554+242 



0601+245 




Fig. 9. — The contour images at several frequen- 
cies for the VLBA and the VERA comparison: 
The images for the four sources are given in each 
column. At the end of each row the array, fre- 
quency, resolution, and session date are listed. 
The resolution for all sources is shown by the cross- 
hatched ellipse for 0556-1-238. The field of view for 
all images is 5.4 x 4.0 mas, and the tick-mark are 
at 0.5-mas intervals. The shows the location 
of the phase center assumed for each source. All 
contour levels are at -5, 5, 10, 20, 50(bold), 75, 
99.5 percent of the peak flux density. 
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/^y=-0.093±0.017 mas/yr 



[] 



=-0.032±0.013 mas/yr 



+ 0547+234 (0) 
* 0554+242 (0) 
□ 0601+245 (0) 



-0.4 



-0.2 0.2 

Time from Session 3 (yr) 



0.4 



Fig. 10.— The Motion of the Cores of 0547+234, 
0554+242 and 0601+245: The position of the eore 
radio components, after removal of the position 
offset, are superposed for the three sources. Both 
43 GHz and 23 GHz points are plotted and the 
source data is indicated by the symbol. Each 
source point is slightly displaced in time in or- 
der to distinguish among the sources. The data 
and error estimates comes directions from the ta- 
bles of the source parameters. The average proper 
motion in X and Y are shown by the straight lines. 
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Fig. 11. — Super-resolution images of 0556+238: 
Images for the five sessions are shown at 43 GHz 
and 23 GHz. The; assumed phase center is show by 
the '+'. A stationary point in this source, derived 
from the apparent core motions of the three target 
sources, is shown by the 'box'. The restoring beam 
at both frequencies, shown by the cross-hatched 
elUpse in the left-most diagram of each row, is 0.2 x 
(iV/S')0.1 mas. The contour levels are at -4, 4, 8, 
16, 32, 50(bold), 80, 99.5 percent of the peak flux 
density. The tick marks are separated by 0.1 mas. 
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Fig. 12. — Frequency-Dependent Structure of 
0556+238: The green contours show the Jan 2010 
8.6-GHz VLBA phase-referencing image. The blue 
contours show the Dec 2009 43 GHz VLBA phase- 
referencing image. The green error bar indicates 
the best estimate of location of the ICRF2 posi- 
tion (see text). The contour levels for both images 
are at 5, 10, 25, 50, 75, 99% of the peak. The tick 
marks are separated by 0.5 mas. 
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Fig. 13. — Frequency-Dependent Structure of 
0547+234: The green contours show the Jan 2010 
8.6-GHz VLBA phase-referencmg uxiage. The blue 
contours show the June 2009 23 GHz VLBA phase- 
referencing image. The green error bar gives the 
estimate of the location of the ICRF2 position for 
this source, derived from the registration of the 
0556-1-238 images. The contour levels for the 8.6 
GHz image are at -1,1,2,4,8,16,32, 50, 75, 99% of 
the peak. The contour levels for the 23 GHz image 
are at -5,5,10,20,40,50,75,99% of the peak. The 
tick marks are separated by 0.5 mas. 
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Fig. 14. — Frequency-Dependent Structure of 
0554+242: The green contours show the Jan 2010 
8.6-GHz VLBA phase-referencing image. The blue 
contours show the Dec 2009 43 GHz VLBA phase- 
referencing image. The green error bar gives the 
estimate of the location of the ICRF2 position for 
this source, derived from the registration of the 
0556-1-238 images. The contour levels for both im- 
ages are at 5, 10, 25, 50, 75, 99% of the peak. The 
tick marks are separated by 0.5 mas. 
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Fig. 15. — Frequency-Dependent Structure of 
0601+245: The yellow contours show the Jan 2010 
8.6-GHz VLBA phase-referencing image. The blue 
contours show the Dec 2009 43 GHz VLBA phase- 
referencing image. The yellow filled ellipse and 
yellow box (size equal to the two-rms error) indi- 
cate two possible locations of the ICRF2 position. 
The contour levels for both images are at 5, 10, 
25, 50, 75, 99% of the peak. The tick marks are 
separated by 0.5 mas. 
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+ 0547 + 234-core 

* 0554+242-core 
□ 0601 +245-core 

★ 0556+238-core 

-J] 



-0.4 



-0.2 0.2 

Time from Session 3 (yr) 



0.4 



Fig. 16. — The Relative Position of the Radio 
Cores: The relative position of the radio cores of 
the four sources, after removal of their constant 
position offset, are superposed for the four radio 
cores. Both 43 GHz and 23 GHz points are plot- 
ted separately, and the source data is indicated by 
the symbol. Each source point is slightly displaced 
in time in order to distinguish among the sources. 
The third session on Jun 8 2009 was used as the 
time origin. 
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Table 1 
Source Parameters 



Name 



RA 



J2000 



Correlated Flxix Density 
DEC 8.6 GHz 23 GHz 43 GHz 

Li S L S L 



0547+234 
0554+242 
0556+238 
0601+245 



05:50:47.39090 
05:57:04.71358 
05:59:32.03313 
06:04:55.12138 



23:26:48.1769 
24:13:55.2986 
23:53:53.9267 
24:29:55.0364 



0.28, 0.25 
0.74, 0.24 
0.56, 0.42 
0.80. 0.34 



0.24, 0.21 
0.49, 0.27 
0.32, 0.30 
0.43, 0.21 



0.18, 0.15 
0.22, 0.09 
0.19, 0.17 
0.10, 0.06 



Note. — ^: S = 200 km baseline; L = 5000 km baseline 



Table 2 
Observing Sessions 



Array Obs. Date Freq (GHz) Resol. e/wxn/s (mas) 



VERA 


18-Apr-2008 


23 


1.0 X 1.0 


VERA 


24-May-2008 


23 


1.0 X 1.0 


VLBA 


20-DCC-2008 


23,43 


0.5 X 1.0, 0.3 X 0.6 


VLBA 


22-Dec-2008 


23,43 


0.5 X 1.0, 0.3 X 0.6 


VERA 


17-Jan-2009 


23 


1.0 X 1.0 


VLBA 


16-Mar-2009 


23,43 


0.5 X 1.0, 0.3 X 0.6 


VERA 


20-Apr-2009 


23 


1.0 X 1.0 


VLBA 


08-Jun-2009 


23,43 


0.5 X 1.0, 0.3 X 0.6 


VLBA 


14-Sep-2009 


23,43 


0.5 X 1.0, 0.3 X 0.6 


VLBA 


12-Dec-2009 


23,43 


0.5 X 1.0, 0.3 X 0.6 


VLBA 


31-Jan-2010 


8.6 


1.5 X 3.0 



Table 3 

Structure Parameters for 0556+238, Component 



Session 




43 


GHz 




23 


GHz 






Intensity [Jy] 


Position [mas] 


Intensity [,Jy] 


Position 


[mas] 




Total 


Peak 


X Y 


Total 


Pealc 


X 


Y 


21-Dec-2008 


0.143 


0.126 ± 0.019 


0.00 ±0.01 -0.00 ±0.02 


0.291 


0.278 ± 0.009 


0.00 ± 0.01 


-0.00 ± 0.02 


16-Mar-2009 


0.178 


0.162 ± 0.013 


-0.00 ± 0.01 0.00 ± 0.02 


0.324 


0.311 ± 0.008 


0.00 ± 0.01 


0.00 ± 0.02 


08-Jun-2009 


0.192 


0.175 ± 0.015 


-0.00 ± 0.01 -0.01 ± 0.02 


0.320 


0.311 ± 0.009 


0.00 ± 0.01 


0.00 ± 0.02 


14-Scp-2009 


0.205 


0.186 ± 0.014 


-0.01 ± 0.01 -0.01 ±0.02 


0.352 


0.340 ± 0.010 


-0.00 ± 0.01 


-0.00 ± 0.02 


12-Dec-2009 


0.217 


0.184 ± 0.017 


-0.01 ±0.01 -0.02 ±0.02 


0.342 


0.330 ± 0.012 


-0.00 ± 0.01 


-0.01 ± 0.02 
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Table 4 

Structure Parameters for 0547+234, Component 



Session 




43 GHz 






23 GHz 






Intensity [Jy] 


Position 


[mas] 


Intensity [Jy] 


Position 


mas] 




Total 


Peak 


X 


Y 


Total 


Peak 


X 


Y 


21-Dec-2008 


0.187 


0.165 ±0.045 


-0.14 ±0.01 


0.24 ±0.03 


0.330 


0.330 ±0.049 


-0.16 ± 0.01 


0.27 ±0.02 


16-Mar-2009 


0.164 


0.153 ±0.027 


-0.16 ±0.01 


0.23 ± 0.02 


0.265 


0.265 ±0.020 


-0.18 ± 0.01 


0.19 ±0.02 


08-Jun-2009 


0.128 


0.115 ± 0.027 


-0.19 ±0.01 


0.21 ± 0.02 


0.199 


0.199 ±0.025 


-0.20 ± 0.01 


0.21 ± 0.02 


14-Scp-2009 


0,162 


0,161 + 0.026 


-0,17 + 0.01 


0,15 + 0,02 


0.227 


0.227 + 0,020 


-0,17 -1- 0,01 


0,13 + 0,02 


12-L)(T-2(I()!) 


0,1!)!) 


1 0,(l,SS 


-0,1,^1 1 0,(12 


0,18 1 


(I.2,S1 


(),2,S1 1 (),1()1 


-0,11 : (),():! 


(1,21 1 0,06 












Table 5 














Structure Parameters for 0554+242 






Session 




43 


GHz 








23 GHz 






Intensity [Jy] 


Position 


[mas] 


Intensity [.ly] 


Position [mas] 




Total 


Peak 


X 


Y 


Total 


Peak 


X 


Y 


Component 


















2I-D0C-2OO8 


0.103 


0.075 ± 0.037 


-0.38 ± 0.02 


-0.36 ± 0.04 


0.133 


0,115 ± 0.025 


-0.42 ± 0.01 


-0,42 ± 0.03 


16-Mar-2009 


0.087 


0.080 ± 0.023 


-0.46 ± 0.01 


-0.39 ± 0.03 


0.148 


0.130 ± 0.019 


-0.44 ± 0.01 


-0.41 ± 0.02 


08-Jun-2009 


0.097 


0.091 ± 0.026 


-0.46 ± 0.01 


-0.40 ± 0.03 


0.151 


0.140 ± 0.020 


-0.46 ± 0.01 


-0.46 ± 0.02 


14-Sep-2009 


0.089 


0.083 ± 0.023 


-0.48 ± 0.01 


-0.45 ± 0.03 


0.146 


0.134 ± 0.019 


-0.46 ± 0.01 


-0.47 ±0.02 


12-Dec-2009 


0.067 


0.062 ±0.026 


-0.47 ±0.02 


-0.45 ± 0.03 


0.126 


0.121 ± 0.041 


-0.45 ± 0.02 


-0.51 ± 0.04 


Component 1 


















21-Dec-2008 


0.113 


0.058 ±0.036 


0.23 ±0.02 


0.14 ± 0.04 


0.348 


0.273 ±0.025 


0.14 ±0.01 


0.08 ±0.02 


16-Mar-2009 


0.117 


0.061 ± 0.022 


0.13 ± 0.01 


0.08 ±0.03 


0.333 


0.249 ±0.018 


0.15 ± 0.01 


0.10 ± 0.02 


08-Jun-2009 


0.124 


0.062 ± 0.025 


0.15 ± 0.02 


0.12 ± 0.03 


0.344 


0.261 ± 0.020 


0.14 ± 0.01 


0.07 ± 0.02 


14-Scp-2009 


0.115 


0.063 ± 0.022 


0.13 ± 0.01 


0.06 ± 0.03 


0.325 


0.245 ± 0.018 


0.15 ± 0.01 


0.06 ± 0.02 


12-Dec-2009 


0.095 


0.050 ± 0.025 


0.14 ±0.02 


0.05 ± 0.04 


0.284 


0.209 ± 0.040 


0.16 ± 0.01 


0.01 ± 0.03 





Table 6 

Structure Parameters for 0601+245 



Session 



43 GHz 



Intensity [Jy] 
Total Peak 



Position [mas] 
X Y 



23 GHz 



Intensity [Jy] 
Total Peak 



Position [mas] 
X Y 



Component 








































2I-D0C-2OO8 


0.053 


0.053 


± 


0.027 


0.84 


± 


0.02 


0.35 ± 


0.04 


0.139 


0.139 


± 


0.031 


0.76 


± 


0.01 


0.25 


± 


0.03 


16-Mar-2009 


0.065 


0.065 


± 


0.024 


0.78 


± 


0.02 


0.32 ± 


0.03 


0.147 


0.147 


± 


0.019 


0.75 


± 


0.01 


0.27 


± 


0.02 


08-Jun-2009 


0.068 


0.068 


± 


0.028 


0.80 


± 


0.02 


0.31 ± 


0.03 


0.127 


0.127 


± 


0.024 


0.78 


± 


0.01 


0.25 


± 


0.03 


14-Sep-2009 


0.058 


0.058 


± 


0.024 


0.78 


± 


0.02 


0.27 ± 


0.03 


0.114 


0.114 


± 


0.024 


0.76 


± 


0.01 


0.25 


± 


0.03 


12-Dec-2009 


0.039 


0.039 


± 


0.033 


0.76 


± 


0.03 


0.21 ± 


0.05 


0.092 


0.092 


± 


0.051 


0.73 


± 


0.03 


0.19 


± 


0.06 


Component 1 








































21-Dec-2008 


0.070 


0.070 


± 


0.027 


0.47 


± 


0.02 


0.24 ± 


0.03 


0.190 


0.190 


± 


0.031 


0.41 


± 


0.01 


0.19 


± 


0.03 


16-Mar-2009 


0.068 


0.068 


± 


0.024 


0.42 


± 


0.01 


0.22 ± 


0.03 


0.165 


0.165 


± 


0.019 


0.37 


± 


0.01 


0.19 


± 


0.02 


08-Jun-2009 


0.061 


0.061 


± 


0.028 


0.44 


± 


0.02 


0.21 ± 


0.03 


0.175 


0.175 


± 


0.024 


0.39 


± 


0.01 


0.17 


± 


0.02 


14-Sep-2009 


0.069 


0.069 


± 


0.024 


0.40 


± 


0.01 


0.19 ± 


0.03 


0.167 


0.167 


± 


0.024 


0.37 


± 


0.01 


0.16 


± 


0.02 


12-DCC-2009 


0.052 


0.052 


± 


0.033 


0.39 


± 


0.02 


0.13 ± 


0.04 


0.156 


0.156 


± 


0.051 


0.34 


± 


0.02 


0.05 


± 


0.04 


Component 2 








































21-Dec-2008 


0.056 


0.029 


± 


0.026 


-0.34 


± 


0.03 


0.13 ± 


0.06 


0.136 


0.080 


± 


0.029 


-0.32 


± 


0.02 


0.10 


± 


0.04 


16-Mar-2009 


0.067 


0.030 


± 


0.023 


-0.44 


± 


0.02 


0.06 ± 


0.05 


0.139 


0.104 


± 


0.018 


-0.44 


± 


0.01 


0.04 


± 


0.03 


08-Jun-2009 


0.064 


0.026 


± 


0.026 


-0.44 


± 


0.03 


0.06 ± 


0.06 


0.125 


0.089 


± 


0.023 


-0.46 


± 


0.02 


0.03 


± 


0.03 


14-Sep-2009 


0.060 


0.022 


± 


0.022 


-0.51 


± 


0.03 


-0.01 ± 


0.06 


0.117 


0.077 


± 


0.023 


-0.50 


± 


0.02 


0.02 


± 


0.04 


12-Dec-2009 


0.047 


0.023 


± 


0.033 


-0.59 


± 


0.04 


-0.06 ± 


0.09 


0.114 


0.061 


± 


0.049 


-0.57 


± 


0.04 


-0.05 


± 


0.08 
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Table 7 

Double Structure Parameters for 0556+238 



Session 






43 GHz 






23 


GHz 






Intensity [Jy] 


Position 


[mas] 


Intensity [Jy] 


Position 


[mas] 




Total 


Peak 


X 


Y 


Total 


Peal; 


X 


Y 


One Component 


















21-Dec-2008 


0.160 


0.125 ± 0.010 


0.00 ±0.02 


0.00 ± 0.03 


0.276 


0.255 ± 0.040 


0.00 ±0.02 


-0.01 ± 0.03 


16-Mar-2009 


0.208 


0.163 ± 0.020 


0.00 ± 0.02 


0.00 ± 0.03 


0.306 


0.280 ± 0.030 


0.00 ± 0.02 


-0.02 ± 0.03 


NE Component 


















08-.Tun-2009 


0.160 


0.160 ± 0.020 


0.01 ± 0.02 


0.03 ± 0.03 


0.280 


0.280 ± 0.040 


0.01 ± 0.02 


0.02 ± 0.03 


14-Scp-2009 


0.180 


0.180 ± 0.020 


0.01 ± 0.02 


0.02 ± 0.03 


0.270 


0.270 ± 0.030 


0.01 ± 0.02 


0.01 ± 0.03 


12-Dec-2009 


0.189 


0.189 ± 0.010 


0.00 ± 0.01 


0.01 ± 0.02 


0.320 


0.320 ± 0.030 


0.01 ± 0.02 


0.02 ± 0.03 


SW Component 


















08-Jun-2009 


0.070 


0.070 ± 0.010 


-0.04 ± 0.02 


-0.11 ± 0.04 


0.045 


0.045 ± 0.015 


-0.04 ± 0.02 


-0.12 ± 0.03 


J. i-Sr|)-2non 


0.018 


n.ni8 1 n.on? 


-n.on 1 0.02 


-n.if) 1 n.03 


n.oon 


n.oGO ■ 0.01.") 


-n.OR 1 0.02 


-0.16 1 n.o.T 


12-DCC-2UOIJ 


U.U7U 


0.U7U 1_ o.uiu 


-O.Ui) ± U.02 


-u.2;i ± o.u:i 


O.IKiO 


0.0(i0 0.01.") 


-0.09 ± 0.02 


-o.iij ± o.o:s 
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